Turtle body size evolution is determined by lineage‐specific specializations rather than global trends

Abstract Organisms display a considerable variety of body sizes and shapes, and macroevolutionary investigations help to understand the evolutionary dynamics behind such variations. Turtles (Testudinata) show great body size disparity, especially when their rich fossil record is accounted for. We explored body size evolution in turtles, testing which factors might influence the observed patterns and evaluating the existence of long‐term directional trends. We constructed the most comprehensive body size dataset for the group to date, tested for correlation with paleotemperature, estimated ancestral body sizes, and performed macroevolutionary model‐fitting analyses. We found no evidence for directional body size evolution, even when using very flexible models, thereby rejecting the occurrence of Cope's rule. We also found no significant effect of paleotemperature on overall through‐time body size patterns. In contrast, we found a significant influence of habitat preference on turtle body size. Freshwater turtles display a rather homogeneous body size distribution through time. In contrast, terrestrial and marine turtles show more pronounced variation, with terrestrial forms being restricted to larger body sizes, up to the origin of testudinids in the Cenozoic, and marine turtles undergoing a reduction in body size disparity after the extinctions of many groups in the mid‐Cenozoic. Our results, therefore, suggest that long‐term, generalized patterns are probably explained by factors specific to certain groups and related at least partly to habitat use.

The evolution of body size is commonly explained in the light of different hypotheses that attempt to elucidate patterns of disparity, such as Cope's rule (or Deperet's rule; Stanley, 1973), a hypothesized tendency lineages to evolve toward larger body sizes (Cope, 1896).
Although directional trends of increasing body size have been identified in a few groups, including medium-to large-bodied mammals (Alroy, 1998) and pterosaurs (Benson et al., 2014), non-directional patterns are present in many other groups Godoy et al., 2019;Laurin, 2004;Moen, 2006). The generalization of hypothesized evolution toward large body sizes has also been questioned, because there is no apparent need for selection to produce ever larger sizes in most lineages (Gould, 1988).
In spite of the presence of hypotheses to explain evolutionary variation of body size in ectothermic vertebrates, few groups have been studied in a comprehensive manner (but see Gearty & Payne, 2020;Godoy et al., 2019;Heim et al., 2015;Smith et al., 2016), with most studies focusing on groups of overall higher metabolic rates, such as mammals and dinosaurs, including birds (Benson et al., 2014Cooper & Purvis, 2010;Cullen et al., 2020;Gearty et al., 2018;Kubo et al., 2019;Raia & Meiri, 2011). Turtles, in particular, have a rich fossil record and relatively stable phylogenetic relations, which provide a reliable framework for macroevolutionary studies. The presence of fossils is particularly important because the inclusion of deep-time data can have a large effect on the outcomes of macroevolutionary analyses.
Considering the relatively low extant diversity of turtles, particularly when compared to mammals, birds, or squamates, the body size disparity of the group is striking. The smallest living testudine, Homopus signatus, has an adult carapace length of about 100 mm, and the largest one, Dermochelys coriacea, reaches more than 2200 mm (Rhodin et al., 2021). Furthermore, fossils display an even broader range of body sizes, including the South American Stupendemys geographicus, with a carapace length of more than 2800 mm (Cadena et al., 2020). This emphasizes the importance of including the available fossil diversity when characterizing patterns of body size evolution in Testudinata (and see Finarelli & Flynn, 2006;Fritz et al., 2013), which has been largely disregarded in most previous attempts (Eastman et al., 2011;Jaffe et al., 2011;Moen, 2006;Uyeda & Harmon, 2014).
Using data from living species only, previous studies presented several hypotheses to explain the observed body size variation of turtles. For example, it has been suggested that such variation is intrinsically related to habitat, with marine species and island tortoises usually possessing larger sizes than freshwater and mainland taxa (Jaffe et al., 2011). This resembles the large body size attained by marine mammals (e.g., Gearty et al., 2018) and, in some aspects, the "island rule" seen in some mammals, reptiles, and birds (Lomolino, 2005). However, Uyeda and Harmon (2014) analyzed turtle body size using unconstrained evolutionary models and suggested that the scenario for turtle optimal body size is more complex than simple differences in habitats, with multiple macroevolutionary body size shifts along the tree. Furthermore, analysis of tortoise (Testudinidae) body size, including fossils, did not find support for the island effect (Vlachos & Rabi, 2018). Moen (2006) tested for evolutionary trends in extant cryptodires, but found no support for directional body size evolution. Moreover, it is not currently clear how temperature influenced long-term body size patterns in turtles, with ongoing discussion about the group following overall rules such as a latitudinal gradient, with larger body sizes seen in colder regions (Angielczyk et al., 2015;Ashton & Feldman, 2003). Yet, except for Vlachos and Rabi (2018) and Sterli et al. (2018), these hypotheses have yet to be tested in a framework including both extinct and extant taxa. In this study, we compiled the largest body size dataset ever assembled for Testudinata, which was used to investigate the tempo and mode of body size evolution in the group, as well as test for possible biotic and abiotic drivers.

| Body size data
Straight-line maximum dorsal carapace length (SCL) was used as a proxy for turtle body size (Jaffe et al., 2011). Aiming to maximize sampling, we also used linear regressions to estimate SCL from the ventral skull length (measured from the rostral tip of the premaxillae to the caudal tip of the occipital condyle) for some specimens lacking carapace. About 7.5% of the SCL data in our dataset was estimated from the ventral skull length. Measurements were collected from photographs (personal archive or the literature), using software ImageJ (Schneider et al., 2012). The final dataset includes body size data for 795 taxa, considerably more than in previous studies (e.g., Angielczyk et al., 2015 = 245 taxa;Jaffe et al., 2011 = 226 taxa;Moen, 2006 = 201 taxa;Vlachos & Rabi, 2018 = 59 taxa). In addition, we also collected habitat preference and chronostratigraphic information for these same taxa using the literature and the Paleobiology Database (PBDB).

| Supertree construction and time calibration
To account for major uncertainties within the phylogenetic relations of the main groups of Testudinata, two informal supertrees were manually assembled using Mesquite version 3.61 (Maddison & Maddison, 2018). These were based on two phylogenetic hypotheses, Evers et al. (2019) and Sterli et al. (2018), hereafter referred to as "Ev19" and "St18," respectively. The most significant differences between the two supertrees are the positions of Protostegidae and Thalassochelyidia (sensu Joyce et al., 2021). Protostegids are stem-Chelonioidea and Thalassochelyidia are stem-Pleurodira in "Ev19," whereas both groups belong to the turtle stem-lineage in "St18" (in which they are originally represented only by Santanachelys gaffneyi and Solnhofia parsoni, respectively). Less inclusive groups were positioned based on several additional hypotheses (Table A1).
Both supertrees include four outgroup taxa (Eunotosaurus africanus, Eorhynchochelys sinensis, Pappochelys rosinae, and Odontochelys semitestacea), which were used for calibration purposes. Each supertree includes a total of 846 taxa, 659 of which are shared with our body size dataset.
Both supertrees were time-scaled using Bayesian inference under a fossilized birth death process (Heath et al., 2014;Stadler, 2010), performed with MrBayes version 3.2.7 (Ronquist et al., 2012). We used R (version 4.0.2; R Core Team, 2021) package paleotree (Bapst, 2012) to create a MrBayes command for timecalibration analyses. The function createMrBayesTipDatingNexus() allows the use of "empty" morphological matrices in clock-less tip-dating analyses (Bapst, 2012;Gearty & Payne, 2020;Godoy et al., 2019). The two supertrees ("Ev19" and "St18") were entered as topological constraints (i.e., for two separate time-scaling analyses) and data on occurrence times (= tip ages) were obtained from the primary literature and supplemented by the PBDB. We used uniform constraints on the tip ages, and the tree age prior was set as a uniform distribution defined between the Kungurian and Roadian stages of the Permian (283.5 and 268.8 million years ago, Ma), as this would represent a maximum possible age for the origin of the group. All other priors were unaltered from the default setting of the createMrBayesTipDatingNexus() function, which were guided by the best practices of Matzke and Wright (2016;see Gearty & Payne, 2020 for more details). Two MCMC runs, with four chains each, were set for 20,000,000 generations, with 25% of the trees discarded as burn-in. Convergence of both runs was verified when values of potential scale reduction factors approached 1.0 and average standard deviation of split frequencies was below 0.01. For both supertrees, we used either the maximum clade credibility (MCC) tree or a set of 10 randomly selected trees from the post-burn-in posterior to perform subsequent analyses.

| Characterizing body size patterns in Testudinata
The entire body size dataset of 795 taxa was used to construct body size through-time plots. Welch's two sample t-tests (Welch, 1947) were used to assess significant changes across different time intervals (i.e., Triassic, Jurassic, Early Cretaceous, Late Cretaceous, Paleogene, Neogene, and Quaternary), focusing on mean body size and disparity, using the standard deviation as a metric of body size disparity. To assess the influence of ecology on the body size distribution, habitat preference information (i.e., terrestrial, freshwater, and marine) was also incorporated into the body size through-time plots. To further test the influence of ecology, we used analysis of variance (ANOVA), performed with R function aov(), as well as the RRPP approach (randomizing residuals in a permutation procedure; , which accounts for phylogenetic dependency (i.e., "phylogenetic ANOVA"), performed with the lm.rrpp() function, from the R package RRPP , and using the MCC tree of each supertree ("Ev19" and "St18").
We also tested for the presence of phylogenetic signal in the body size data using the R function phyloSignal() (Keck et al., 2016), using 10 randomly selected trees from the posterior distribution of trees of both supertrees. We used 1000 replicates and estimated To further characterize body size evolution within Testudinata, we used maximum likelihood to estimate ancestral body sizes under Brownian motion (BM), using the fastAnc() function of the R package phytools (Revell, 2012). Inferred ancestral sates were performed with both the complete supertrees (i.e., using the MCCT trees with all 659 taxa, including fossils and extant species) and a subtree with only extant taxa (i.e., dataset reduced to 312 taxa).

| Testing for the presence of Cope's rule
To test if Cope's rule played an important role in turtle body size evolution, we fitted different evolutionary models to our body size data in both supertrees. To account for temporal and phylogenetic uncertainties, 10 time-scaled versions of each alternative supertree ("Ev19" and "St18") were used.
We fitted four uniform phenotypic models to our data, starting with the uniform BM model, in which body size undergoes an unconstrained, single-rate random walk along phylogenetic lineages, resulting in diffusive evolutionary expansion (Felsenstein, 1973(Felsenstein, , 1985Freckleton & Harvey, 2006). This pattern is consistent with several possible causes, including genetic drift or wandering adaptive optima (Felsenstein, 2003), between which genetic drift seems a less likely explanation at macroevolutionary scales. The model has two parameters: sigma squared (σ 2 ), which indicates evolutionary rate, and the root state of the trait at time zero, sometimes represented by X(0) (Felsenstein, 1973).
We also fitted three other uniform models: (1) the "mean_trend" (or "drift") model, which is a modification of the BM model that incorporates a parameter (μ) describing an uniform directional trend along all branches of the phylogeny (Pagel, 2002); (2) the EB model (also known as "ACDC model"; Blomberg et al., 2003), in which lineages experience a burst of rapid increase in trait variation in the beginning of their evolutionary history, followed by a deceleration (Harmon et al., 2010); (3) the Ornstein-Uhlenbeck (OU) model, which incorporates attraction of trait values (represented by the α parameter) toward an optimum (θ) (Butler & King, 2004;Hansen, 1997). In the case of the OU model fitted here, the parameters (α and θ) were not allowed to vary along the tree.
We also fitted 13 non-uniform trend-like models to our data.
Unlike uniform "mean_trend" model, these multi-regime models allow the μ parameter-the amount of directional change in a trait through time (Hunt & Carrano, 2010;Pagel, 2002)-to vary along the tree in temporal or node shifts ("time-shift model" and "node-shift model," respectively). We fitted "time-shift" models (which allow shifts in all branches after a determined point in time) allowing the number of shifts to vary from one to three; and "node-shift" models (which allow shifts in some branches) allowing the number of shifts to vary from 1 to 10, resulting in a total of 13 multi-trend models (3 "time-shift" and 10 "node-shift" models). In this study, both uniform trend and multi-trend models were fitted as a representation of the Cope's rule, given that it is described as a multi-lineage directional trend toward larger sizes (Cope, 1887(Cope, , 1896Stanley, 1973).
Moreover, among mammals, the foundational example of Cope's rule, directional body size evolution is only present in some lineages (e.g., Alroy, 1999), more consistent with a "node-shift" model, with multiple independent origins of directional evolution.

| Influence of paleotemperature
We used regressions to test for the possible influence of global paleotemperature on turtle body size (795 taxa). As a proxy for paleotemperature, we compiled δ 18 O data (lower δ 18 O values indicate higher environmental temperature) from two different sources. First, we used tropical isotopic data collected in tropical regions by Prokoph et al. (2008), who assembled isotopic information from marine organisms, extending from Precambrian to recent. Furthermore, we also used global paleotemperature data from Zachos et al. (2008), which compiled information about isotopic ratios in foraminifer shells from the Maastrichtian to the recent. We tested for correlation between both temperature curves and our body size indices, including maximum, minimum, and mean body size, as well as body size disparity (= standard deviation).
Correlation between body size data and paleotemperature was initially assessed using ordinary least squares. In addition, to avoid potential issues created by temporal autocorrelation, we used generalized least squares with a first-order autoregressive model incorporated (tsGLS; Fox & Weisberg, 2018), using the R package nlme (Pinheiro et al., 2022). The data were divided into time intervals, using approximately equal-length (~9 million years) stratigraphic time bins (from Mannion et al., 2015). For each time bin, we calculated body size indices (disparity [= standard deviation], maximum, minimum, and mean body size) and weighted mean δ 18 O values using R package disparity (Guillerme, 2018).
The figures presented by this study were made using the R packages palaeoverse (Jones et al., 2023) and deeptime (Gearty, 2023).

| Inferred ancestral states
Inferred ancestral states based on both supertrees ("Ev19" or "St18";  Vlachos and Rabi (2018), which included 59 taxa (23 living and 36 extinct). Moreover, our ancestral body size estimate for Testudines (345 mm) is smaller than that estimated by Sterli et al. (2018; 359 mm), whereas our Testudinata ancestral body size is much smaller (396 mm, in comparison to 570 mm estimated by Sterli et al., 2018). More importantly, however, some similar patterns are observed between our results and those by Sterli et al. (2018), such F I G U R E 1 Ancestral body sizes (log 10 maximum dorsal carapace length in millimeters) mapped onto Testudinata phylogeny ("St18"), with (a) complete tree and (b) extant-only subtree. Ancestral body size for different taxonomic groups (c); small gray dots indicate all taxa within that lineage; colored triangles represent the ancestral estimates of stem-groups; circles represent ancestral estimates of the crown groups; and squares represent ancestral estimates of the crown groups without the fossil taxa. Gray area indicates sizes between 200 and 500 mm. The numbers indicate the same groups in a-c.
as the body size increase seen in Podocnemididae and the decrease seen in Kinosternidae.

| Model fitting
The AICc scores for all the evolutionary models fitted to the turtle trees and body size data show an overwhelmingly stronger support (i.e., lower AICc values) for the uniform OU model, even when compared to the non-uniform multi-trend models. Consistently, this stronger support for the OU model was found when using both "Ev19" and "St18" topologies (Table 1). These results rule out the presence of trend-like processes (either uniform or multi-trend) in the body size evolution of Testudinata, at least when the entire tree is considered.

| Correlation with paleotemperature
No significant correlation was observed between global paleotemperature and mean, minimum, or maximum body size of Testudinata through time (Table A2). However, we did find a significant, but weak correlation between body size disparity (= standard deviation of body sizes) and paleotemperature (Table 2; Figure 2), with disparity increasing at lower temperatures (i.e., higher δ 18 O values). Although significant, this correlation is relatively weak when palaeotemperature data from Prokoph et al. (2008) are used, but it becomes slightly stronger when using data from Zachos et al. (2008), which is restricted to the Late Cretaceous-Recent time interval. This may suggest that the influence of environmental temperatures on turtle body size was stronger during the Cenozoic.

| Differences among ecological habitats
Freshwater is the most common habitat occupied by turtles from the Jurassic onwards, and this category is therefore the most influential to the aggregated pattern of turtle body size variation through time. No significant changes in either disparity or mean body size of freshwater turtles occurred since the Late Cretaceous (Figure 3b,d).
A significant increase in mean body size in freshwater turtles was identified between the Early (234 mm) and the Late Cretaceous to 443 mm; p-value = .02865). However, body size disparity is low, as expected from the low number of species, and the group is mostly represented by medium-or large-bodied forms (Figure 3a).
The first marine turtles appeared at the end of the Jurassic, with body sizes like those of other groups (Figure 3c). There is a noteworthy (although not significant; p-value = .3872) increase in the mean body size of marine turtles from the Jurassic to the Early Cretaceous, TA B L E 1 Results of model-fitting analyses, depicting model parameters and AICc scores for the models fitted to our body size dataset of Testudinata (log 10 maximum dorsal carapace length) and 10 time-calibrated trees for each of the two initial supertree topologies ("St18," based on the hypothesis of Sterli et al., 2018, and   and a small drop in the Late Cretaceous (Figure 3c,d), but the latter interval witnessed the highest variability in body sizes among marine turtles in the series (Figure 3c). After the K-Pg transition, the range of body sizes decreased substantially, although the mean remains similar. After that, from the Neogene to the Quaternary, a significant increase in mean body size (p-value = .003958), from 413 to 1036 mm, was detected for marine turtles (Figure 3c,d).
We found that ecological habitat is significantly linked to body size in turtles using ANOVA (Table 3). However, when the phylogenetic structure of the data is taken into account (using phylogenetic ANOVA; Table 3), this association is only significant for the freshwater and terrestrial ecological categories. Moreover, for all trees tested here, we found a strong phylogenetic signal in body size data (λ > 0.856; p < .001).

| Body size patterns and ancestral estimates: The effect of including fossils
Despite previous controversies (e.g., Patterson, 1981), it has become increasingly clear that the paleontological record is crucial and Figure A3). The latter has been affected for most lineages assessed here, but to different degrees (Figure 1), with no directional effect on the estimates (Figure 1). This differs from the pattern seen in mammals, for which ancestors have considerably larger body sizes when estimated using fossils (e.g., Finarelli & Flynn, 2006;Finarelli & Goswami, 2013). For mammals, these results were explained by the widespread occurrence of directional evolution toward large body sizes from small ancestors (Alroy, 1999;Smith et al., 2010), a pattern not observed in turtles (see below).

| Cope's rule and directional trends of body size evolution
We found no evidence for directional patterns of body size evolution in turtles (Table 1), given that none of our trend-like models (either the uniform trend or the multi-trend models) received TA B L E 2 Results of time series generalized least squares (tsGLS) and ordinary least squares (OLS) regressions using turtle body size disparity and paleotemperature data. Coef. more support than the uniform OU model. This result is consistent with previous investigation by Moen (2006), which also found no support for Cope's rule when analyzing extant cryptodires (even though directional evolution may be difficult to detect on extantonly datasets; Finarelli & Flynn, 2006;Schnitzler et al., 2017;Slater et al., 2012). Therefore, our results provide no support for the hypothesis that Cope's rule explain the evolution of large body sizes in Testudinata. They also add to the growing evidence that directional body size evolution is rare among vertebrates Godoy et al., 2019;Huttenlocker, 2014;Sookias et al., 2012), with the exception of mammals and pterosaurs (Alroy, 1998;Benson et al., 2014), once again challenging the generality of Cope's rule.

| Influence of environmental temperature on turtle body size evolution
The relationship between abiotic factors and body size has been extensively studied, with distinct vertebrate groups being differently affected by them, especially when it comes to comparing endothermic and ectothermic organisms (Angielczyk et al., 2015;Angilletta et al., 2004;Ashton & Feldman, 2003;Mousseau, 1997;Partridge & Coyne, 1997;van der Have & de Jong, 1996;Van Voorhies, 1996).
Large-scale trends such as Bergmann's rule (i.e., the tendency of having larger body sizes at higher latitudes within a species) may play an important role for within-species variation of body size in endotherms (Ashton et al., 2000;James, 1970;Zink & Remsen, 1986), and may explain patterns of maximum size during mammal evolution (Saarinen et al., 2014). Yet, results for ectothermic reptiles are less consistent (Angielczyk et al., 2015;Ashton & Feldman, 2003;Mousseau, 1997).
We evaluated the correlation between turtle body size distributions and paleotemperature variation through time. In general, no significant influence of temperature on mean, minimum, or maximum body size in the group was found (Table A2). Similar results were reported for crocodylomorphs, for which no significant correlation between paleotemperature and body size (mean, maximum, and minimum values) was found when the entire group is analyzed, even though a strong association between both variables is observed when only the crown group is considered (Godoy et al., 2019). Therefore, although our results indicate no overall influence of paleotemperature on the through-time distribution of turtle mean, maximum, and minimum body sizes, we cannot rule out an influence of environmental temperatures at smaller temporal and phylogenetic scales. Indeed, environmental temperature has been a commonly proposed explanation for body size variation in different turtle clades and species, particularly affecting disparity, diversity, or distribution of less inclusive groups (e.g., Böhme, 2003;Ferreira et al., 2018;Georgalis & Kear, 2013;Vitek, 2012).

F I G U R E 2 (a)
Through-time patterns of Testudinata body size disparity (standard deviation of log 10 maximum dorsal carapace length in millimeters) and paleotemperature (δ 18 O isotopic data from Zachos et al., 2008) during the last ~70 Ma. Error bars were calculated by bootstrapping the disparity data 500 times. δ 18 O is used as proxy for paleotemperature and is inversely proportional to temperature. (b) Linear regression (OLS) between turtle body size disparity and δ 18 O data (regression results shown in Table 2). We did find a significant correlation between paleotemperature and turtle body size disparity (= standard deviation) during the Cenozoic, with periods of higher disparity associated with lower temperatures (or higher δ 18 O values; Figure 2). We suggest three potential explanations for this seemingly counterintuitive result.
First, low temperatures might have restricted niche availability for turtles and, consequently, driven body size specialization-toward larger or smaller body sizes-to avoid competition, as seen in some extant lineages (Cunha et al., 2020;Pritchard, 2001). Conversely, colder and dryer environments could have increased availability of coastal habitats (by sea level drops), which has been associated with higher diversification rates in turtles (Thomson et al., 2021).
Higher species richness might have also led to higher disparity levels in body size. Finally, the significant correlation might be an arti-

| Body size evolution and ecological habitats
Our ANOVA results (Table 3)  such habitats (Jaffe et al., 2011), with different species, closely related or not, inhabiting several disparate freshwater environments (Bonin et al., 2006). For instance, different pleurodiran and cryptodiran lineages acquired resistance to estuarine or brackish water (Agha et al., 2018;Bower et al., 2016). Also, closely related taxa occupying the same areas are known to avoid competition through body size divergence, such as extant podocnemidids (Cunha et al., 2020) and trionychids (Pritchard, 2001).
Terrestrial and marine turtles, on the other hand, are represented by fewer lineages, with more restrict evolutionary histories. The earliest turtles were terrestrial, ranging from medium-to large-sized during the Mesozoic (Figure 3a-d). Meiolaniformes is the only of these stem-lineage groups to survive until recently (until the Holocene; , displaying large to gigantic body sizes, especially after the Mesozoic. Testudinids-the only extant lineage of terrestrial turtles-appeared in the fossil record during the Paleogene and remained relatively small until at least the end of that period (Figure 3c). The Eocene-Oligocene witnessed a peak of diversity, related to the origin of crown-group Testudinidae (Lourenço et al., 2012;Vlachos & Rabi, 2018), after which the group spread from Eurasia to most of the world. Therefore, the recent high variation in body size within tortoises might be also related to the expansion of occupied habitats (as in freshwater turtles), also associated with specific diversification and extinction dynamics within the group (e.g., Joos et al., 2022).
Extant marine turtles (Chelonioidea) exhibit low disparity, but remarkably large body sizes (Figure 3b), which might be related to morphological adaptations to a pelagic lifestyle, given that other groups of sea turtles (e.g., Thalassochelydia and Bothremydidae) are not as strongly associated with larger sizes and were probably not pelagic.
The large size of marine turtles might also be explained by either physiological constraints (e.g., thermoregulation; Mrosovsky, 1980) or the need for higher dispersal abilities associated with migration TA B L E 3 Results of ANOVA and PhyloANOVA exploring the relationship between Testudinata body size (log 10 maximum dorsal carapace length in millimeters) and different ecological categories: terrestrial (n = 145), freshwater (n = 563), and marine (n = 87 Note: Pairwise comparisons between ecological categories also shown. Abbreviations: df, degrees of freedom; MS, mean squares; SS, sum of squares. *Significant at alpha = .05. (Jaffe et al., 2011). It has been previously proposed that thermoregulation and other physiological aspects (e.g., lung capacity while diving) play an important role in determining the larger body sizes of aquatic mammals and reptiles Davis, 2014;Gearty et al., 2018;Gearty & Payne, 2020;Gutarra et al., 2022;Pyenson & Vermeij, 2016;Williams, 2001), by posing a minimum body size limit on these species. However, in the case of marine turtles, the fossil record shows that smaller species also existed in the past (Figure 3c,d), with Santanachelys gaffneyi from the Early Cretaceous of Brazil as one of the oldest and smallest sea turtles (200 mm;Hirayama, 1998).
Therefore, in the case of Testudinata, perhaps the lower body size limit imposed by physiological constraints was not as strict as those inferred for other secondarily aquatic tetrapods (e.g., mammals). On the other hand, the shell and the necessity to lay eggs on land possibly pose constraints on the maximum body sizes achieved by marine turtles , which are, in general, smaller than other Mesozoic marine reptiles and extant cetaceans Smith & Lyons, 2011).
Finally, intrinsic factors might also influence body size evolution in turtles. Sterli et al. (2018), for example, suggested that a reduction in size in Mesochelydia (the clade including all post-Triassic turtles)which is confirmed by our ancestral state estimates (Figure 1)-could be explained by paedomorphic processes, which are also evidenced by other morphological traits. The shell might also constrain the maximum body size of turtles inhabiting terrestrial and semiaquatic environments (Golubović et al., 2017;Lyson et al., 2014). It could hamper turtles from attaining sizes as large as giant mammals and dinosaurs due to a different relation between weight and body size.
Moreover, minimum body sizes in turtles are overall larger than that of the smallest lissamphibians, squamates, mammals, and birds.
Endothermy could explain the smaller sizes of mammals and birds (Lovegrove, 2017), but not of lissamphibians or squamates. Hence, it is possible that the shell also imposes a lower body size limit to testudinatans. We did not find support for a general trend-like process leading to larger body sizes, discarding Cope's rule as an explanation for body size evolution in turtles. Also, we did not find a significant influence of paleotemperature on mean, maximum, and minimum body size. Although we found a significant, moderate correlation between temperature and body size disparity through time, this association might be an artifact caused by a join constant increasing of disparity and continuous drop in temperatures during the Cenozoic.

CO N FLI C T O F I NTE R E S T S TATE M E NT
The authors have no conflicts of interest to declare.

O PE N R E S E A RCH BA D G E S
This article has earned an Open Data badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://osf.io/a45hv/ ?view_ only=ffa79 065ef f54a8 48f17 168c1 1f3e4f6.

DATA AVA I L A B I L I T Y S TAT E M E N T
The complete dataset, supertrees, and scripts used for the analy-

F I G U R E A 2
The "St18" time-calibrated supertree, represented by the 50% majority rule tree. Blue bars represent the 95% highest posterior density (HPD) age ranges for each node.

F I G U R E A 3
Ancestral body sizes (log 10 maximum dorsal carapace length in millimeters) mapped onto Ev19 supertree, with (a) complete tree and (b) extant-only subtree. Ancestral body size for different taxonomic groups (c); small gray dots indicate all taxa within that lineage; colored triangles represent the ancestral estimates of stem-groups; circles represent ancestral estimates of the crown groups; and squares represent ancestral estimates of the crown groups without the fossil taxa. Gray area indicates sizes between 200 and 500 mm. The numbers indicate the same groups in a-c. Coef. Abbreviations: Coef., coefficient; interc., intercept; p.R 2 , Nagelkerke pseudo R-squared; SD, standard deviation; Var., variable.